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Abstract — We describe new methods for deciding the stability 
of switching systems. The methods build on two ideas previously 
appeared in the literature: the polytope norm iterative construc- 
tion, and the lifting procedure. Moreover, the combination of 
these two ideas allows us to introduce a pruning algorithm 
which can importantly reduce the computational burden. 

We prove several appealing theoretical properties of our 
methods like a flniteness computational result which extends 
a known result for unlifted sets of matrices, and provide 
numerical examples of their good behaviour. 

I. Introduction 

A switching discrete time linear dynamical system in M" 
is denned by 



x t +\ 



(1) 



where ^# = {Ai} ieJ! ( J? is a set of indices) is a set of 

n x n real matrices that describes the system, and xq is the 

initial state. So, a trajectory is uniquely defined by an initial 

condition xq and a switching signal o{t ) : N — > , which 

represents the sequence of matrices ruling the dynamics. 

Switching systems constitute an important family of hybrid 

systems. They have been the subject of a great research 

effort in recent years (see [10], [15], [19], [28] for general 

introductions and applications in systems and control). Of 

particular importance in many applications is the worst-case 

stability of these systems. It is well known that their stability 

is ruled by the joint spectral radius of the set of matrices 

^f, which we now define: 

For each k = 1,2,..., consider the set ,^ k {^() of all 

possible products of length k whose factors are elements 

of that is @> k {J?) = {Ai l ...A ik | i\,...,i k € J] and 

set &(^) = |J & k (y$() t° be the multiplicative semigroup 
k>\ 

associated with 

Definition 1 (Joint Spectral Radius - jsr [27]): If || ■ || is 
any matrix norm on 
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The joint spectral radius (jsr) of ^ is defined as 

p{J() = lim p k (J(). (2) 

The joint spectral radius does not depend on the manix 
norm chosen thanks to the equivalence between matrix norms 
in finite dimensional spaces. It is a natural generalization of 
the classical notion of spectral radius of a matrix to a set 
of matrices. It is also known to be very hard to compute 
[29]. Yet, it is of high importance in practice, because it 
characterizes the stability of switching systems: 

Theorem 1: (e.g. [15]) The system (HJ converges to the 
origin for any initial point xq and any switching signal a if 
and only if 

p{Jt) < 1. 

One could come up with another natural generalization of 
the spectral radius to a set of matrices: 

Definition 2 (Generalized Spectral Radius - gsr [7]): 
Let p(-) denote the spectral radius of an n x n-matrix, 
consider p~ k (jX) := sup Pe ^^ p(P) x l k , k G N and 
define the generalized spectral radius (gsr) of as 

p(^#) = limsup ~p k (^). (3) 

In the case of bounded sets ^ Daubechies and Lagarias 
conjectured the equality of gsr and jsr, which was then 
proven by Berger and Wang [2]. From now on we shall 
restrict ourselves to finite sets of matrices, and we call the 
joint spectral radius of a finite set ^# the quantity 

pM) = pM) = pM). 
It is not difficult to see that for any bounded set of matrices, 

Vk,p p~(^) k < p{J() < Ppi-d). (4) 

The joint spectral radius has found a great number of 
applications in several fields of Engineering, Mathematics, 
and Computer Science. See [15] for diverse applications. 
Recently, several methods from very different flavors have 
been proposed to provide approximations of this quantity 
[1], [9], [20], [28]. In particular, there has been a growing 
interest in the study of the so-called extremal norms (also 
called Barabanov norms) [5], [17], [21]. We now detail this 
idea. 

A. Extremal norms 

An important characterization of the joint spectral radius 
p(^#) of a matrix set is the generalization of the follow- 
ing well-known formula, p(A) = inf [| . || eOp |K||> where Op 
denotes the set of all possible operator norms. 



In order to state this characterization, we define the norm 

of the set ^# = {A,} ((S ,y as 

||^||= ft (>f)= max ||A,-1|. 

Proposition 1 (see e.g. [15]): The joint spectral radius of 
a bounded set ^ of n x n-matrices is given by 

pL#) = inf |L#||. (5) 

|| - 1| eOp 

We remark that the right-hand side of (01 is often referred 
to as the common spectral radius of the set ^ . 

Given a set ', it is important to establish if the infimum 
in © is indeed a minimum. We have the following definition. 

Definition 3 (Extremal norm): A norm || • ||* satisfying 
the condition 

\\Jt\\,=p{Jt) (6) 
is said to be extremal for the set 

For sets of matrices with jsr equal to one, a sufficient 
condition for the existence of an extremal norm is that the 
product semigroup generated by ^ is bounded and the set 
is irreducible. 

An extremal norm is a norm whose unit ball is invariant 
under the action of the matrices in the set. This confers 
a practical interest to the analysis of such objects: if one 
manages to generate a norm whose unit ball is an invariant 
set, then this directly allows to bound the jsr: p < 1. 

Remark 1: We consider the particular question p < 1 
without loss of generality because one can divide all the 
matrices by a constant r and ask whether the newly obtained 
set has joint spectral radius smaller than one. By homogene- 
ity of the Definition [2] the minimal r such that the answer 
is yes provides the actual value of the joint spectral radius. 
This is the strategy that we will follow in this paper. 
Methods based on the generation of an invariant polytope 
have been proposed in [12]— [14], [23]. A great advantage of 
these methods is that they can stop in finite time, and provide 
an exact computation of the jsr of the set of matrices. Indeed, 
if, besides the knowledge of an invariant polytope, one knows 
a product whose spectral radius is equal to one, then (by Eq. 
@) one has a proof that the exact value of the jsr is one. 

B. Invariant cones 

Along an other line of research, methods have been 
proposed, which base their analysis on the existence of a 
common invariant cone for the set of matrices, or which first 
transform the set of matrices so that the new matrices share 
an invariant cone ( [3], [16], [22], [25], [26]). In this paper, by 
cone we mean a convex, closed, pointed, and nondegenerate 
cone. Let K C R" be a cone with the apex at the origin (see 
e.g. [4]). Any such cone defines a partial order in R": we 
write x y (x >k y) for x — y& K(x — y& intK). 

The cone K is an invariant cone for the matrix A if AK C K. 
In this case we say that A is nonnegative and write A >k 0. 
If K is invariant for all matrices of some set . //. then it is 



said to be an invariant cone for that set. We now present a 
simple algebraic operation, that enables one to construct a 
new set of matrices without essentially changing the JSR. 
Moreover it turns out that for any set of matrices the 
newly constructed set has an invariant cone: the set S^" of 
semidefinite positive matrices. 

Definition 4: Given a vector v 6 C", its semidefinite lifting 
is defined as 

v = Re(vv*). 

Also, the semidefinite lifting of a matrix is the linear operator: 

A:,9"l^,9"l, Av=Av. 
We will note a < b (resp. a>^b) for a <y« b (resp. a >yn b). 

The following important result relates the spectral radii of 
a set of matrices and the corresponding lifted set [25]. 

Theorem 2: The following equality holds true for any 
bounded set of matrices 

p(^)=p{Jtf. 

Thus, if one wants to compute the joint spectral radius^of 
he can as well compute it for the transformed set ^ . 
Our goal is to exploit the invariant cone property of the 
lifted set in order to compute the joint spectral radius by 
determining an extremal norm for M . As we will see, it 
will be sufficient to find an invariant set inside the invariant 
cone. A convenient way of describing such a set is given in 
the following definition. It is reminiscent of the convex hull 
operation, but, besides being defined by a set of points, it 
also depends on the geometry of an underlying cone. 

Definition 5 (Conical convex hull and Conitope): Given 
a cone K and a (possibly infinite) set of vectors U C K, its 
conical convex hull is defined as 

Coni({/) = {xeK:3ye Conv(C7) : x < K y}. 

We say that a set S C K is a conitope if there is a finite 
set U ~ {uj}j = i such that 

S = Com(U) = {xeK :3A,i,...,^,Ai>0,^A,- = l: 

x^<X\ii\A vXiU(\. (7) 

Conitopes can in turn define a norm on the corresponding 
cone in the following way. 

Lemma 1: Let K be a cone. Let U be a set of points 
in K, such that Conv(t/) H intK ^ 0. Then, for any x G K, 
\x\ v = inf{A > : x/X e Coni(t/)} is a valid norm on K. 
We call such norms conitope norms. 

We present a last definition which characterizes minimal sets 
of vertices defining a conitope. 

Definition 6 (Essential system): We say that U C K is an 
essential system of vertices if for any U' dU,U' =/= U, 

Coni(C/')^Coni(C7). 

II. Plan 

In this work, we show how to combine the two above 
mentioned approaches: our goal is to find conitope norms 
which are extremal for a given set of matrices. We will show 
that this combination allows to derive innovative efficient 
algorithms for computing the joint spectral radius. As a 



byproduct, the same ideas also provide an approximation 
algorithm for the joint spectral subradius of a set of matrices. 
However, we focus here on the joint spectral radius computa- 
tion for the sake of clarity. In the next section, we present the 
algorithms that unify the two approaches. Then, in Section 
IIVI we prove that the algorithms are correct and converge to 
the required value. In Section [V] we exhibit a large family of 
sets of matrices for which our algorithms provably converge 
in finite time. In section [VT] we give evidence for the good 
performance of our method on numerical examples. In the 
conclusion, we point out open questions and generalizations 
that naturally come up from our methods. 

III. The algorithms 

A. Computing the joint spectral radius 

We start by giving the following useful definition. 

Definition 7 ( Spectrum-maximizing product): If ^# is a 
bounded set of complex n x n-matrices, any matrix A £ 
& t (JZ) satisfying 

P (A) 1 / t =p t (^)= P {j£). (8) 

for some t > 1 is called a spectrum-maximizing product 
(SMP) for J(. 

The procedure of Algorithm 1 can be summarized in three 
steps, which are analogous to those used in the classical 
polytope algorithms ( [12]— [14]). The first step is to look for 
a candidate spectrum maximizing product A £ £P t (^) for 
some t > 1 and to scale the set as 

Ji' = J(jp(A) x l\ 

which satisfies p (./#') > 1. 

The secondj3tep is to look for an invariant convex set for 
the lifted set i.e., the unit ball of a conic extremal norm. 

If the procedure succeeds, then we can conclude that 
p(^#') = 1, that is, the jsr of our initial set of matrices is 
actually equal to p(A) l / t . 

Some escaping criteria to detect that A is not an SMP, 
can be used in analogy to those used for classical polytope 
algorithms (see e.g. [6], [11]). 

Algorithm 2 builds an extremal conitope in a dynamic 
manner, without reinitializing the algorithm, and pruning the 
products to consider thanks to the partial order implied by 
the cone y". 

B. Two algorithms 

We present here our two algorithms. They both generate 
an upper bound Y and a lower bound C which converge 
towards the joint spectral radius. 

We denote as L.E.(A) the leading eigenvector of a matrix 
A, implicitely assuming this is unique. In the case it is not 
unique we will specify this. 

Algorithm 1: the conitope method 

Algorithm 1: 

. Preprocessing: 

• Find a product A of length t such that p (A) 1 1' is large (if 
needed, perform an exhaustive search for increasingly 
long products) 



. C = p(A) 1 /' 
. vo:=L.E.(A) 
. Uq := vb 
. J(' := JtjC 

- while int 5„ + nspan([/o) = 

* Let Uq be an essential system of vertices of 
Com(U U^'U ) 

- U:=U 

• Main loop 

- Repeat until B = 1 

* W := jf'U 

* B :=max{|w| £/ : w eW},Y = B-C 

* Let U be an essential system of vertices of 
Coni([/UW) 

* D := max{|w|j/ : w £ Uq} 

* If D < 1 

• A is not an SMP: go back to preprocessing 

- Return Coni(£/) as an invariant conitope. 

• End 

Algorithm 2: the dynamical procedure 

We now present a modification of the algorithm which 
is allowed by the fact that the matrices share a common 
invariant cone. The idea is that the cone induces a partial 
order relation between points, which allows us to disgard 
some of them, if they are dominated by another point. 
This algorithm is a dynamical version of the previous one: 
while trying to prove that a certain product A is an SMP, at 
the same time it tries and finds another better product A'. 
Once it finds such a product, the algorithm simply continues 
its computation, except that the stored points in V are scaled 
as if the algorithm had started with A' . 

Algorithm 2: 

. Preprocessing: 

. C = max {p (A) : A £ 

. Jt = Je)c 

. U :=I 

. U:=U 

• Main loop 

- Repeat until B = 1 

* W := JW 

* B := max {\w\u:weW}, Y=BC 

* Let U be an essential system of vertices of 
Coni(t/UW) 

* L = max A {p(A) 1 /'W : A£/(z),z£W}(/(z) is 
the set of subproducts of the product that maps 
Uq on z. 1(A) is the length of the product A.) 

* if L > 1 

■ For all ueU : u:=u/(C 1 ^) (l(u) is the length 
of the product A such that AUq = u.) 

- Return Coni([/) as an invariant conitope. 
Remark 2: The algorithms differ by two facts: first, in 

the initialization part, the first one starts with a leading 



eigenvector of a candidate SMP, and the second one with the 
identity matrix; second (and mainly), in the first algorithm, 
one reinitializes the conitope, by restarting from the new 
leading eigenvector, while the second algorithm keeps all 
the previous vertices and continues with them. These two 
modifications are independent from each other, but the main 
idea is that in the first algorithm, starting with a leading 
eigenvector enables us to stop in finite time (indeed, it is 
known that the leading eigenvectors of the permutations of 
the SMP are all on the boundary of any extremal norm). 
For the second algorithm, we start with the identity matrix 
because we do not really hope to converge in finite time, 
but we are more interested in obtaining rapidly converging 
bounds than by a finite time termination. 

Finally, in algorithm 1 as presented here, we do not search 
for a better spectral radius among all the subproducts in 
f(z) for all z G W, like we do in algorithm 2. However, 
one could perform this search as well in Algorithm 1, which 
would allow to potentially find a better SMP without having 
to resort to an exhaustive search when coming back to the 
preprocessing phase. 

IV. Bounds and convergence results 

In this section, we show that our methods are correct, 
which means that they provide valid lower and upper bounds 
on the joint spectral radius, and that moreover, these bounds 
converge towards the actual value asymptotically. 

Proposition 2: Let Jt be a set of matrices that share an 
invariant cone K. Let U be a set of points in K, such that 
Conv([/) DintK ^ 0. Then, we have the following valid upper 
bound on the joint spectral radius: 

p(^)<pu(^), (9) 

where 

Pu{Jt) = sup sup {|A«|f/}. (10) 

Proof: By Equation^yJ), for any e > 0, therejydsts a 
t and a product A G £?{Jjf) such that p(A) > (p(JC) - ej . 
Since the matrices in Jt share an invariant cone K, by the 
generalized Perron-Frobeniy« Theorem, there exists a vector 
u G K such that Au = (p(Jt) — e)'u. Let us take the vector 
u such that \u\jj = 1. 

Now, Au G pu^JtJV (as implied by Eq. ( flOl i) together with 
Au= (p — e)'u imply 

Pu{~^) > P -£• 

since this must hold for any e we obtain ©. ■ 
Since algorithm 1 starts with a set of vectors t/o, which 
are rank one matrices, the set U actually contains only rank 
one matrices during the whole run of the algorithm. Hence, 
it is not clear that one obtains for sure a valid norm (as 
required to go out of the while loop in the algorithm). The 
next proposition shows how to deal with this problem: 

Proposition 3: If the matrices in J( do not have a com- 
mon invariant real subspace, then Algorithm 1 generates in 
finite time a set U such that Conv(C/) nintS"" ^ 0. 
Moreover, if the algorithm does not generate such a set U 



after n iterations, then (with the notations of Algorithm 1) 
S = span ^ a (^#)(Re(vo)) is an invariant real subspace, and 
one can project the matrices on S and on 5 X and iterate the 
algorithm. 

Proof: By contraposition, suppose that 

Conv{(^(^)v )} n mt,y n+ = 0. This implies that 
there exists a real vector x such that 

VA G &>{^C),x*A(V )A*x = 0. 

Then, letting 

Vo = Vi+!V 2 , V; G K", 

we have that 

VA G ^(^),JC*A(viVi+v 2 v|)A*jc = 0, (11) 
which implies that 

x*Av i v* i A*x = 0, i = 1,2 

and &(Jt)vi is a non trivial invariant subspace for Jt ' . 

Moreover, Equation ([Til l holds if and only if it holds for 
all A G ^(Jf), k = 1,2, . ..n. Thus, if it still holds at step 
n, one can project Jt on (for instance) &(Jt)v\ and on its 
orthogonal space, and reiterate the algorithm on the two sets 
of smaller matrices. ■ 

In the proof of our main theorem below, we also need the 
following technical proposition: 

Proposition 4: [8] Let || • || be a matrix norm in C 
induced by the vector norm • | . There is an absolute constant 
c(n) > such that for all z G C, \z\ = 1, and all A G 
C" , | |A 1 1 < 1 , there is an eigenvalue A of A such that 

1 1 A, | <c\Az~z\ 1/n . 
Theorem 3: Both algorithms 1 and 2 are exact, i.e., both 
the upper bound Y = C * B and the lower bound C asymp- 
totically converge towards the joint spectral radius. 
The proof can be found in the extended version of this paper 


Remark 3: If the set Jt is not irreducible, then we can 
still use Algorithm 2, but we cannot converge for sure 
towards an invariant conitope, because such a conitope might 
not exist. However, the lower bound is still valid, and still 
converges towards to JSR. For the upper bound, one has to 
resort to the classical estimate p < sup^g^dlAjl ''''). 

V. FlNITENESS RESULTS 

A. Preliminary results: asymptotic simplicity and extremal 
polytope norms 

We denote here the normalized set 




Now we report the existence results proved in [12], which 
involve the normalized trajectories 

Sr\Jt',x\ = {x} U {Px I P G .<?{Jt')} 

1 Lifted polytope methods for the asymptotic analysis of matrix semi- 
groups,¥veprmU 2012, by the same authors. 



where = \J S? k {JC). First we observe that all the 

k>o 

cyclic permutations of a product P have the same eigenvalues 
with the same multiplicities. Thus, if P = A;, . . .A,, is an 
SMP for a set M ' , then each of its cyclic permutations 

A (j . . . Ajj Ai k , . . . A, t+1 , S = 1, . . . , K 1, 

still is an SMP for Jf . 

Definition 8 (Asymptotically simple set): A nondefective 
bounded set ./# of complex n x n-matrices is called asymptot- 
ically simple if it has a minimal SMP P with only one leading 
eigenvector (modulo scalar nonzero factors) such that the set 
X of the leading eigenvectors of ^ is equal to the set of the 
leading eigenvectors of P and of its cyclic permutations. 

We remark that the asympotic simplicity of means that 
the set X of the leading eigenvectors of ^ is -cyclic, i.e., 
for any pair (x,y) e X x X there exist a, j3 G C with 

\a\-\p\ = l 

and two normalized products P',Q' £ such that 

y = aP'x and x = j3 Q'y. 

Theorem 4: [12] Assume that a finite set . // of complex 
n x «-matrices is nondefective and asymptotically simple. 
Moreover, let x ^ be a leading eigenvector of and 
assume that span! SF\M' ,x\ ) = C". Then the set 



(12) 

is finite (modulo scalar factors of unitary modulus). As 
a consequence, there exist a finite number of normalized 
products P[, . . . ,P' S G g?{JC) such that 

S?\Jt\x\=&szo{\x,P[x y ...,P' s x\\, (13) 

so that ,x\ is an invariant balanced complex polytope. 

B. Asymptotic simplicity and extremal conitope norms 

We assert that if a set of matrices is asymptotically simple, 
not only it admits an extremal polytope norm as claimed by 
Theorem[4] but its lifted set also admits an extremal conitope 
norm. The proofs are to be found in the extended version of 
this paper. We start with a technical lemma. 

Lemma 2: If a sequence of vectors s, £ C" converges to 
zero, the set Coni ({.?,•}) is a conitope. 

Theorem 5: If a finite set of matrices ^ is_Jrreducible 
and asymptotically simple, then its lifted set j$ admits an 
extremal conitope. Moreover, there exists such an extremal 
conitope whose all vertices are rank one matrices. 

Remark 4: 

Note that if one applies the polytope algorithm in [12], [14] 
and then lifts the vertices {v;} and considers the correspond- 
ing conitope, one does not obtain in general an invariant set 
for the lifted set of matrices. This is a key point in order 
to understand that the two algorithms are fundamentally 
different (see the following Example [3] in Section [VTb . 



VI. Numerical examples 

We provide here some illustrative examples demonstrating 
the effectiveness of the methods proposed in this article. 
Example 1: Consider the set of matrices: 
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-1 
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-i 
-i 
-i 



-i 
-i 
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This set of matrices has a joint spectral radius 
p(^T) =p(A 2 ) « 1.779. 

The classical BCP algorithm terminates in 4 steps and 
finds a BCP which is the convex hull of 16 vertices. 
Algorithm 1 terminates after 2 steps and finds an invariant 
conitope with only 7 vertices. 

Example 2: Our next example presents a set of complex 
matrices. 




In this case p(jK) = p(AiAiA 2 AiA 2 ) ~ 2.2401. 

The classical BCP algorithm requires 20 steps to find 
a BCP with 65 essential vertices. Algorithm 1 terminates 
after 8 steps and the invariant conitope has only 10 essential 
vertices. 

Example 3: The following set of matrices has an invariant 
polytope convjv,}, but it appears that simply lifting its 
vertices does not provide an invariant conitope: Coni({v;}) 
is not invariant for . ^ = {Ai,A2} with 



A, = 




A? = 




The spectrum maximizing product is P=A\A 2 and has a real 
leading eigenvector. An extremal real centrally symmetric 
polytope has been computed by the algorithm described in 
[13] and gives 6 vertices: 

v\=x, v 2 =A' 1 vi, v 3 =A 2 vi, 

v 4 =A' 2 v 2 , v 5 =A 2 v 3 , v 6 =A' 1 v 5 , 

where A\ = A,/p(P), i = 1,2. 

The simple algebraic test that w = A' 2 v 2 does not belong to 
the conitope shows that computing the set of lifted vertices 
yields a conitope which is not invariant. 

VII. Conclusion 

In this paper we introduced a new method for deciding the 
stability of switching systems. By combining the lifting tech- 
nique and the extremal polytope technique, we have shown 
that one can in some cases outperform the other algorithms 
in the literature. We have also provided a sufficient condition 



for the existence of an extremal norm for the lifted matrices 
("an extremal conitope"). 

A huge advantage of our algorithms is that they allow to 
disgard in a systematic way many products in the semigroup, 
thanks to the partial order relation induced by the semidef- 
inite cone. This can greatly reduce the computation time in 
some situations. Another advantage is that the algorithms 
can terminate in finite time, and provide a proof that a 
particular product is actually an SMP, by exhibiting a convex 
set (the extremal conitope) which is a valid extremal norm. 
The same ideas apply for the joint spectral subradius, and 
provide the first efficient stopping criterion that allows to 
compute exactly the joint spectral subradius for general sets 
of matrices (i.e. with negative and positive entries). We delay 
for further work the empirical study of these methods for the 
computation of the joint spectral subradius. 
Also, we leave for further research the question of whether 
similar methods could be applied in order to approximate 
other joint spectral characteristics, like the so-called p- 
radius, or the Lyapunov exponent. Indeed, it has been shown 
recently that other methods based on semidefinite program- 
ming are useful for computing these quantities as well (e.g. 
[18], [24]). 

The main computational burden, at every step, is in the 
computation of the quantity \Av\y, i.e., the semidefinite 
program allowing to compute the conitope norm of the new 
vertices. If the matrices do not share a priori a common 
invariant cone, one must work in the lifted space in order to 
have an invariant cone (that is, ,5 fi n +), but then one has to 
resort to semidefinite programming to represent the conitope 
on a computer. We leave it as an open question whether it 
is possible to apply similar ideas to general matrices with a 
simpler cone, on which the computations are less expensive. 
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Appendix 

Proof of Theorem [3] 
Proof: In the following, p denotes p(^#) = p{Jt) 2 . 

Lower bound With the notations of Algorithm 1, let us 
call the scaled set of matrices: 

J(' = JK/C. 

We want to prove that the lower bound C converges towards 
the actual value of the Joint Spectral Radius. Thus, for 
Algorithm 1, it suffices to show that if p («/#') > 1 (i.e. if A is 
not an SMP), the main loop terminates in finite time, and the 
algorithm will hence look for longer products by returning 
to the preprocessing phase; since we have the well known 
alternative definition (01 for the JSR, the lower bound will 
converge towards p. 

From Proposition [3] we can suppose that we have n vectors 
V\,...,v n £ ^(^#')vo linearly independent. 

Since p(^K') > 1, it is possible to find arbitrary large 
vectors satisfying this property. In particular, they can be 
such that 

£(v})/rt h v . 

Hence, the algorithm will terminate, because the last 
equation implies that Uq C Coni(f/) at some time, which 
is exactly the termination condition in Algorithm 1 . 

For algorithm 2, recall that for any t and any norm | • |, 

p{Jff < sup{|A«| : \u\ = ljei"}. 

Thus, there exists a vector u S Uq, an infinite prod- 
uct . . . A t . . . A i , and a sequence f,-, , f; 2 , . . . such that 
\A t . ...A\u\ > p'' . Let us now consider the normalized 
sequence u' t = u t /\u t \, u t = A t . . .A\u. This is an infinite 
sequence on the boundary of the unit ball of the norm. By 
compactness of this boundary, there exists a subsequence s,- 
such that the u' s . converge to a vector y : \y\ = 1. So for any 
e, there must exist Sj such that for all s, > sj, 

\u{. — u'\ < e. 



Setting z = u' s .,Bi =A Sj . . .A Sj+ i(\u Sj \/\u Si \), we get matrices 
Bj e 8?{J() such that 

\B iZ -z\ < e, 

and we can conclude by Proposition @] that Bt has an 
eigenvalue X such that |1 — A| = C^e 1 /"), which implies that 

p{A Si ...A s . +l ) l ^-*^ P . 

As a conclusion, Algorithm 2, which checks the spectral 
radius of all the products of the type A tj ...A tj , will find 
candidates SMP whose averaged spectral radius tends to p. 

Upper bound Again from Proposition [3] we can suppose 
that ./# has no real invariant subspace. This implies [15, 
Theorem 2.1] that the semigroup & (^(^ / p)j is bounded. 

Therefore, denoting by C t the lower bound C at step f, 
we have U, = {Jt '/C t )*U tends to t/„ = &>{.Jt/p)Uo in 
Hausdorff distance. 



Hence the quantity sup 



\Au\u t tends towards 



sup IAmIj/^ = p(,/, 



Appendix 



Proof of Lemma |2] 
Proof: 

Define W = span^Si}^)), k ~ dim(W), and take a set 
{ri,...,r^} C {s;};>o of full rank. Then, all the s, can be 
written as Y,j=i,...,k^j r j- 

Take now a vector Sj with i large enough such that 



{r h ...,r k } C {si,...,,y,-_i}, 



and such that, denoting Sj = ^j r ji < l/^ 2 - 

y'=l 

We claim that 



^d£(or;)A. 

7=1 

Indeed, taking an arbitrary vector igR", we have 



x*s i s* i x = &hh2){ x * r h r *h x ) 
ii J2=i 

h J2=l 

< k 2 {\/k 2 fmwt\r*x\ 2 



(14) 



(15) 



(£(r^)/*)x. 



Since it is true for any ieR", (TT~5T > follows by definition. 

As usual we indicate as F= Re(ss*). Equations (fl4l i and 
(fT5T > together imply that 

Sj eConi(si,...,^_i), 

and then, 

Coni({Si,...,i5}) =Coni({ii,...,^_i}). 
Now, take i such that 

k 

VZ >i,3Ai,...,A jt : < l/k 2 ,s, = ^kjrj. 

Thus, for i large enough we have 

Coni({si,...,S;}) =Coni({Si,..., 

which implies that the limit of this set actually is a conitope, 
since the limit is reached for a finite number of vertices. ■ 



Appendix 

Proof of Theorem [5] 
Proof: We assume without loss of generality that p (^#) = 
p(~#) = 1. Let us denote X the finite set of leading 
eigenvectors (modulo scalar nonzero factors) of Take a 
particular v EX and denote V the set of trajectories starting 
at v = Re(vv*) under arbitrary products of matrices in M : 

V = {A il ...A h ^):A ij e^ Vj}. 

We now define 

S = Coni(V). (16) 

Since ^# is irreducible, again by [15, Theorem_2.1], the 
set V is bounded, and S is an invariant set for ^ . The last 
follows by the fact that if U is an invariant set, then Coni(E/) 
is also an invariant set. 

We will prove that there exist a finite number of vectors 
Vi,...,V; G V such that S = Coni({vi,. . -,vi}), and so S is 
an invariant conitope for . 

Suppose the contrary. 
Then, there must be an infinite set of vectors V' CV such 
that 

Wi£V',VigCom(V'\vi). (17) 

Moreover, by compactness of the set S we can suppose that 
V' is a subset of a single trajectory of v, that is, there exist 
indices i'i , Z2 5 - - - such that the set 

V' = {A it ...A h v:t >0} 

contains an infinite set satisfying Equation ilTi . (For a 
detailed proof of this fact see for instance [12, Lemma 5.19].) 

We note v, = A,, ...A (1 v, and correspondingly 
Vt = A( t . . .A,jV. Since this set is infinite, there must 
be an index ?o such that v, E X, v, +i $X. 

Now, by compactness of Coni(y'), there must exist a 
subsequence of v, that converges to a vector w. This vector w 
is different from zero because, by Equation (fTTI i and Lemma 
|2] the sequence A,, . . . A (1 v is bounded away from zero. Then, 
since 

v t = Re(v(V*), 

this must be the case also for the original (unlifted) matrices: 
hence 

3t < h < t 2 < ■ ■ .tk • • • : A iti . . • A,-, o+2 v, 0+1 w, 

and then w EX. Now, by asymptotic simplicity, there exists 
a matrix P G such that Pw = v tQ . This implies 

that PAi tk ■ ■ -Ait + 2 V <o+! — * v 'o' an( ^ nence ' premultiplying the 
above equation by A,- ( +] , we obtain that 

A i, 0+ i P \ ■ ■ - A k 0+ 2 v t +l -> y ro+l> 
and v, 0+ i is in X, which gives a contradiction. ■ 



